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Gear Crack Propagation Path Studies - Guidelines for Ultra-Safe Design 

David G. Lewicki 

U.S. Army Research Laboratory 
National Aeronautics and Space Administration 
Glenn Research Center 
Cleveland, Ohio 
david.g.lewicki@grc.nasa.gov 

Design guidelines have been established to prevent catastrophic rim fracture failure modes when considering gear tooth 
bending fatigue. Analysis was performed using the finite element method with principles of linear elastic fracture mechanics. 
Crack propagation paths were predicted for a variety of gear tooth and rim configurations. The effects of rim and web 
thicknesses, initial crack locations, and gear tooth geometry factors such as diametral pitch, number of teeth, pitch radius, and 
tooth pressure angle were considered. Design maps of tooth/rim fracture modes including effects of gear geometry, applied 
load, crack size, and material properties were developed. The occurrence of rim fractures significantly increased as the 
backup ratio (rim thickness divided by tooth height) decreased. The occurrence of rim fractures also increased as the initial 
crack location was moved down the root of the tooth. Increased rim and web compliance increased the occurrence of rim 
fractures. For gears with constant pitch radii, coarser-pitch teeth increased the occurrence of tooth fractures over rim 
fractures. Also, 25° pressure angle teeth had an increased occurrence of tooth fractures over rim fractures when compared to 
20° pressure angle teeth. For gears with constant number of teeth or gears with constant diametral pitch, varying size had 
little or no effect on crack propagation paths. 


Introduction 

Effective gear designs balance strength, durability, 
reliability, size, weight, and cost. However, unexpected gear 
failures may occur even with adequate gear tooth design 
(Ref. 1). In order to design an extremely safe system, the 
designer must ask and address the question “what happens 
when a failure occurs.” With regard to gear tooth bending 
fatigue, tooth or rim fractures may occur. For aircraft, a 
crack which propagates through a rim would be catastrophic, 
leading to disengagement of a rotor or propeller, loss of an 
aircraft, and possible fatalities (Refs. 2 and 3). This failure 
mode should be avoided. A crack which propagates through 
a tooth itself may or may not be catastrophic, depending on 
the design and operating conditions. Also, early warning of 
this failure mode may be possible due to advances in modem 
diagnostic systems (Ref. 4). 

Fracture mechanics has developed into a useful 
discipline for predicting strength and life of cracked 
structures. Linear elastic fracture mechanics applied to gear 
teeth has become increasingly popular. Among the earliest, 
fracture mechanics was applied to gear teeth to simulate 
crack propagation, compute threshold loads, estimate stress 
intensity factors, and calculate tooth life (Refs. 5-7). 
Researchers at Tohoku University in Japan performed a 
series of analyses and experiments to determine the effect of 
residual stress on crack initiation and propagation (Refs. 8 
and 9). In addition, a comprehensive, self-contained analysis 
package to refine the spur gear bending fatigue theory using 
fracture mechanics was developed (Ref. 10). 

The stress intensity factors are the key parameters to 
estimate the characteristics of a crack. Analytical methods 
using weight function techniques to estimate gear tooth stress 
intensity factors have been developed (Refs. 11 and 12). 
Numerical techniques such as the finite element method and 


boundary element method have also been studied (Refs. 13 
and 14). Based on stress intensity factors, fatigue crack 
growth and gear life predictions have been investigated 
(Refs. 15-18). In addition, gear crack trajectory predictions 
have been addressed in a few’ studies (Refs. 19-25). 

The objective of the current study is to develop design 
guidelines to prevent catastrophic rim fracture failure modes 
when considering gear tooth bending fatigue. Analysis was 
performed using the finite element method with principles of 
linear elastic fracture mechanics. Crack propagation paths 
were predicted for a variety of gear tooth and rim 
configurations. The effects of rim and web thicknesses, 
initial crack locations, and gear tooth geometry factors such 
as diametral pitch, number of teeth, pitch radius, and tooth 
pressure angle were considered. Crack trajectories are 
presented for the variety of cases studied along with design 
maps indicating gear tooth or gear rim fracture modes. It 
should be noted that the current study investigates the 
likelihood of tooth or rim fracture assuming an initial crack 
is present. The absolute probability of fracture should 
include crack initiation, but is beyond the scope of this wxirk. 

Gear Modeling 

Basic gear tooth geometry data was input to a tooth 
coordinate generation computer program. The tooth 
coordinate generator program used the method of Ref. 26 to 
determine the tooth coordinates. The output was tooth 
coordinate and rim coordinate data which defined a single- 
tooth sector of a gear. This output was used by a 
commercially available pre- and post-processing finite 
element analysis software package (Ref. 27). This package 
created the finite element mesh of the complete gear. The 
mesh was then imported to the FRANC (FRacture ANalysis 
Code) computer program. 
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FRANC is a general purpose finite element code for the 
static analysis of cracked structures (Ref. 28). The program 
is designed for two-dimensional problems, uses principles of 
linear elastic fracture mechanics, and is capable of analyzing 
plane strain, plane stress, or axi-symmetric problems. Eight- 
node quadrilateral or six-node triangular elements can be 
used. Among the variety of capabilities, a unique feature of 
the program is the ability to model a crack in a structure. The 
program uses a method called "delete and fill" to accomplish 
this. To illustrate, consider a finite element mesh of an 
uncracked structure. The user would first define an initial 
crack by identifying the node of the crack mouth and 



Fig, 1 . Crack modeling scheme using finite element 
method, (a) user-defined initial crack, (b) final mesh 
of initial crack, (c) predicted crack propagation path. 


coordinates of the crack tip (Fig. la). FRANC would then 
delete the elements in the vicinity of the crack tip, insert a 
rosette of quarter-point, six-node triangular elements around 
the crack tip to model the inverse square -root stress 
singularity, then fill the remaining area between the rosette 
and original mesh with conventional six-node triangular 
elements (Fig. lb). The user would then run the finite 
element equation solver to determine nodal displacements, 
forces, stresses, and strains. Mode I and mode II stress 
intensity factors, K t and K Ih respectively, can be calculated 
using a variety of methods. (As a refresher, mode I loading 
refers to loads applied normal to the crack plane which tend 
to open the crack. Mode II refers to in-plane shear loading.) 
The stress intensity factors quantify the state of stress in the 
region near the crack tip. In the program, the stress intensity 
factors can also be used to predict the crack propagation 
trajectory angles, again using a variety of methods. 

A further unique feature of FRANC is the automatic 
crack propagation capability. After an initial crack is 
inserted in a mesh, the program simulates crack propagation 
as a number of straight line segments. For each segment (or 


step), the program solves the finite element equations, 
calculates the stress intensity factors, and calculates the 
crack propagation angle. The program then places the new 
crack tip at the calculated angle and at a user-defined crack 
increment length. The model is then re-meshed using the 
"delete and fill" method described above. The procedure is 
repeated a number of times as specified by the user. Fig. lc 
shows the predicted crack propagation path of a gear tooth. 
In this example, the predicted crack trajectory was after 29 
steps, i.e., the crack trajectory was approximated by 29 line 
segments. In previous studies, gear crack propagation paths 
calculated from the FRANC computer program were 
validated from experimental tests (Refs. 19, 24). Such a 



Fig. 2. Comparison of predicted gear tooth crack 
propagation paths with experimental results, 
P=predicted, E=experiments (Ref. 24). (a) m b - 3.3, 
(b) m b = 1.0, (c) m b - 0.5. 


validation is shown in Fig. 2. Here, the effect of rim 
thickness (expressed as backup ratio, m b , defined in the next 
paragraph) on crack path was explored. Notches were placed 
in the fillet region of the test gear teeth and run in a fatigue 
test rig until tooth or rim fracture occurred. The FRANC 
program was also used in these studies to model the gears. 
Initial cracks were inserted in the tooth fillets (corresponding 
to the notch locations of the test gears) and propagated as 
described above. As seen from the figure, the program was 
successful in predicting the crack paths. 

A typical finite element gear model used in the current 
study is shown in Fig. 3. The mesh shown is for an 
uncracked gear. The gear design is the baseline used in the 
current study. The design parameters are: 28 teeth, 
8 diametral pitch, 1.75-in pitch radius, 20° pressure angle, 
and a 0.25-in face width. The model had 2255 plane stress, 
8-node, quadrilateral elements and 7122 nodes. For 
improved accuracy, the mesh was refined in the upper 
portion of the model (this is the region where cracks will be 
inserted). The tooth load was placed at the highest point of 
single tooth contact (HPSTC), normal to the surface. 
Although the tooth load changes in magnitude and direction 
in actual gear operations, a static analysis with the load at 
the HPSTC has given accurate results with respect to crack 
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Fig. 3. Typical finite element model of an uncracked 

gear; 28 teeth, 8 diametral pitch, 1 .75-in pitch radius, 

20° pressure angle, m b - 1.0. 

propagation analysis (Ref. 24). Four hub nodes at the gear 
inner diameter were fixed to ground for boundary 
conditions. The material used was steel. In addition, slots 
were incorporated in the model to model thin-rim gears. The 
model shown has a backup ratio, m b =l.O. The backup ratio is 
defined as the rim thickness, b, divided by the tooth height, h 
(Fig. 4a). As stated before, crack propagation angles are 
determined from the calculated stress intensity factors. In the 
current study, the stress intensity factors were determined 
from the finite element method nodal displacements and 
forces using the J-integral method (Ref. 29). In addition, the 
crack propagation angles were determined from the stress 
intensity factors using the maximum tangential stress theory 
(Ref. 30). 

Effects of Backup Ratio and Initial Crack Location on 
Crack Propagation 

Gear models with backup ratios m b of 0.5, 0.6, 0.7, 0.8, 
0.9, 1.0, 1.1, 1.2, and 1.3 were investigated. These models 
were based on that shown in Fig. 3, but with various slot 
heights to give the appropriate rim thicknesses. Also, the 
effect of initial crack location, 6 0 , was investigated. The 
location of the initial crack is defined in Fig. 4b. 0q defines 
the location of the initial crack mouth on the tooth fillet or 
root region with respect to the pitch radius. 

The effect of initial crack location on crack propagation 
path is shown in Fig. 5 for a backup ratio of m b - 1.0. Initial 
crack lengths of 0.010 in were individually inserted and 
propagated in the models until the cracks reached either 
tooth or rim boundaries. This took from 25 to 49 steps, 
depending on the case, using a crack increment length of 
0.010 in. For Fig. 5, the initial crack location angles varied 
from 00 = 6 8 to 120 °. 00 = 88 ° corresponded to the root 
centerline between the tooth teeth (index 5 on Fig. 5). 
00 = 104 ° corresponded to the location of the largest tensile 
stress for an uncracked gear of this design (index 8 on 
Fig. 5). Note that for 00=68 to 78°, rim fractures occurred. 
For 00=83 to 120°, tooth fractures occurred. The mode I and 




Fig. 4. Definition of terms, (a) backup ratio, m b -b!h , 
(b) initial crack location angle, 6 0 . 


mode II stress intensity factors as a function of crack length 
for the m b - 1.0 gear crack models of Fig. 5 are shown in 
Fig. 6. These are for an applied tooth load of Q - 364 lb. For 
clarity, not all the initial crack location cases are shown in 
the figure. For the rim fracture cases (00=68 to 78°), only the 
00=68 case is shown for mode I. The other rim fracture 
cases, however, had nearly identical responses. Also note 
that the K r magnitudes for these cases were rather low. For 
the tooth fracture cases, the magnitude of Ki increased as the 
initial crack location, 0 0 , increased. Also, the K[ magnitudes 
had a significant increase toward the later portion of the 
propagation. Here, the teeth were nearly fractured off. The 
mode II stress intensity factor for 00=83 is shown in Fig. 6b. 
Although not exactly the same, K u for the other initial crack 
cases was in the -1 to 3 ksiVin range. The significant 
observation here is that the mode I stress intensity factors are 
much greater in magnitude than the mode II. This implies 
that the crack propagation paths are smooth, continuous, 
and, in most cases, rather straight with only a slight 
curvature. This matches that seen in field experience for gear 



Fig. 5. Effect of initial crack location on crack 
propagation path, m b = 1.0, 1) 0 y =68°, 2) 0 O =73°, 
3) 00=78°, 4) 00=83°, 5) 0 O =88°, 6) 0 O =94°, 

7) 00=99°, 8) 00=104°, 9) 0 o =lO9°, 10) 0 O =1 14°, 
11 ) 00 = 120 °. 
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Crack Length, in 

Fig. 6. Stress intensity factors from gear crack 
propagation studies, m b = 1.0, 2=364 lb. (a) mode I, 
(b) mode II. 


tooth bending fatigue. (Again, as a refresher, the magnitude 
of K i determines the crack propagation rate while K u 
determines the crack propagation direction. For A7/= 0, the 
crack propagates in a straight path. For K If ^0 and much 
smaller than K h the crack propagates in a slightly curved 
path. The crack propagation angle is a function of K u / AT/). 

Fig. 7 shows the effect of backup ratio and initial crack 
location on crack propagation. For the baseline gear design 
parameters and a backup ratio of m b - 1.3, tooth fractures 
occurred for all initial crack locations (Fig. 7i). As the 
backup ratio decreased, the transition from tooth fractures to 
rim fractures occurred at a larger initial crack location angle. 
Thus, for thinner rims, the occurrence of a rim fracture 
significantly increased. Fig. 8 shows the mode I stress 
intensity factors for a variety of backup ratios. Fig. 8 a is for 
an initial crack location of 0^104° and Fig. 8 b for 0^88°. 
Again, these corresponded to the location of the largest 
tensile stress for an uncracked gear and the location of the 
root centerline, respectively. For 0^=1 04°, the K ( magnitudes 
were nearly identical for backup ratios m^=0.8 to 1.3, For 
thinner rims, the K { magnitudes decreased as the backup 
ratio decreased. This was due to the increased compliance in 
the gear rims, which reduced the tensile stress in the tooth 
fillet region. Rim fracture occurred for m b = 0.5 and tooth 
fractures occurred for m b > 0.6 for the 0^=1 04° cases. For 
0o=88 0 and for crack lengths less than 0.2 in, the K I 
magnitudes were slightly less than those for the 0o= 104°. 
This occurred because the crack tip regions were located 
further down the fillet toward the root where the tensile 
stress was lower. An exception to this was the ^=0.5 case 
where the Kj magnitudes were about the same for 0 o =88 and 




Fig. 7. Effect of backup ratio and initial crack location 
on propagation path. 

104°. For crack lengths greater than 0.2 in, larger crack 
lengths were required for the 0 ^= 88 ° cases compared to the 
00=104° cases to reach final fracture. Rim fracture occurred 
for m^0.7 and tooth fractures occurred for 0.8 for the 
00 = 88 ° cases. 

Fig. 9 summarizes the crack fracture mode in a design 
map. Here, the crack failure mode is plotted as a function of 
both backup ratio and initial crack location. Note that the 
design space is divided into three regions: 1 ) tooth fractures, 
2) rim fractures, and 3) no crack propagation due to 
compression at the initial crack tip. Again, for the baseline 
gear design, 00=104° is the location of the largest tensile 
stress in the tooth fillet. That would be the probable location 
of crack initiation as long as no material defects are 
considered. Thus, one would want to have a backup ratio 
m}> 0.6 to have only tooth failures at that location. On the 
other hand, if the potential for crack initiation at other 
locations is considered, a backup ratio of m^l.3 should be 
used to minimize the probability of rim failures. 

It should be noted from Figs. 6 and 8 that the mode I 
stress intensity factors at the start of crack propagation are 
rather low for lower values of 0 O . Considering this, the 
design map of Fig. 9 can be adjusted using the stress 
intensity factor threshold concept. The stress intensity factor 
threshold, AK th , is largest value of the mode I stress intensity 
factor in which no crack propagation would occur. It is a 
material property and can be derived through standard 
fatigue crack growth tests. Table I shows the normalized 
mode I stress intensity factors as a function of backup ratio 
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Mode I stress intensity Mode I stress intensity 

factor, Kj (ksi\/in) factor, Kj (ksiVin) 



Crack Length, in 

Fig. 8. Stress intensity factors from gear crack 
propagation studies, £>=364 lb (a) 6 0 = 104° 
(maximum fillet stress), (b) 6> 0 =88° (root 
centerline). 
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Fig. 9. Effect of backup ratio and initial crack 
location on crack failure modes, T=tooth fracture, 
R=rim fracture, C=compression (no crack 
propagation). 
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location on crack failure modes, £>=500 lb, 
a =0.030 in, AK rh =5 ksiVin; T=tooth fracture, 
R=rim fracture, N=no crack propagation. 


Table I. Normalized mode I stress intensity factors, Kj IQ (in~ 3/2 ). 


Backup 
ratio, m B 

120 

114 

109 

Initial crack location, 0 0 (deg) 
104 99 94 88 83 

78 

73 

68 

63 

1.3 

14.58 

16.21 

16.68 

15.99 

13.94 

11.28 

8.86 

6.92 

5.21 

3.47 

1.91 


1.2 

14.54 

16.12 

16.55 

15.85 

13.80 

11.17 

8.79 

6.91 

5.24 

3.54 

1.99 


1.1 

14.47 

16.01 

16.40 

15.69 

13.65 

11.06 

8.75 

6.95 

5.34 

3.66 

2.11 


1.0 

14.38 

15.86 

16.22 

15.50 

13.48 

10.97 

8.75 

7.05 

5.51 

3.87 

2.29 


0.9 






10.93 

8.88 

7.29 

5.82 



1.19 

0.8 





13.20 
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1.48 
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13.91 

15.18 
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14.79 

13.22 

11.38 

9.82 

8.52 

7.16 

5.38 

3.45 

1.81 

0.6 



15.17 

14.76 

13.60 

12.28 







0.5 

13.30 

14.49 

15.00 

15.11 

14.84 

14.29 

13.33 

12.22 

10.62 

8.14 

5.30 

2,98 
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and initial crack location for an initial crack of a=0.030 in. 
The normalized stress intensity factors were derived by 
dividing the stress intensity factors from the finite element 
analysis, K h by the applied load, Q. Since linear elastic 
fracture mechanics is assumed, one can scale the normalized 
stress intensity factors for any value of applied load, then 
compare the results to the stress intensity factor threshold. 
Fig. 10 is a modified design map as an example for an 
applied load of 0=500 lb, a stress intensity factor threshold 
AK Jh -5 ksiVin (this is a typical value for AISI 9310 steel, the 
current standard material in aerospace drive system 
applications), and a crack size of a=0.030 in. For many of 
the cases, the mode I stress intensity factors were less than 
the stress intensity factor threshold, and thus, no crack 
propagation occurred. For the conditions of Fig. 10, a 
backup ratio of 7?t^>0.8 should be used to ensure no rim 
failures will occur. This approach of using the stress 
intensity factor threshold concept is probably the most 
realistic, since cracks initiating at low 6 0 conditions are 
rather rare in field experience. However, the design map 
becomes more complex since it is dependent on gear 
geometry, applied load, crack size, and material properties. 

Effect of Fillet Geometry on Crack Propagation 

Fig. 1 1 shows the same basic gear tooth shape with two 
different fillet designs. The first is a standard fillet 
(Fig. 11a). The second is an increased fillet (Fig. 11b), 
which was used as the baseline design in the previous 
section. The standard fillet was derived by increasing the 
number of teeth of the cutting tool in the gear tooth 
generation process (Ref. 26). Fig. 12 shows the design map 
for the effect of tooth fillet on crack propagation. The 
increased fillet slightly increased the proportion of tooth 
fractures over rim fractures. Although not shown, the 
increased fillet had an additional benefit of reducing tensile 
stress in an uncracked gear. 

Effect of Rim/Web Compliance on Crack Propagation 

To first investigate rim compliance effects, a partial 
finite element model of the baseline design was developed 
(Fig. 13). The model was a four-tooth partial model of the 
baseline design (28 teeth, 8 diametral pitch, 1,75-in pitch 
radius, 20° pressure angle) for a backup ratio of m b = 0.9. The 
standard fillet design was used (Fig. 11a). The edge of the 
rim as well as the inner radius were fixed to ground for 
boundary conditions. Although m b =0.9 is not considered a 
thick-rimmed gear, the boundary conditions used tended to 
make the rim extremely non-compliant. The tooth load was 
placed at the HPSTC, normal to the surface. Fig. 14 
compares the crack propagation paths of the baseline slotted 
gear (Fig. 3, but for m b =0.9 and with the standard fillet) to 
the partial-model gear of Fig. 13. The partial model gear had 
tooth fractures for both 6o=S 3 and 88° while the slotted gear 
had a rim fracture for 00=88°. The conclusion reached was 
that increased rim compliance (such as with thin-rimmed 
gears) leads to more rim fractures. 

To investigate web compliance effects, a full, non- 
slotted finite element model of the baseline design was 
developed (Fig. 15). The model had 28 teeth, 8 diametral 



Fig. 11. Gear tooth shapes, 8 pitch, 28 teeth, pitch 
radius 7^=1. 7 5 in. (a) standard fillet, (b) increased 
fillet. 



Fig. 12. Effect of tooth fillet on crack failure modes. 



Fig. 13. Finite element model for rim compliance 
effect on crack propagation; 28 teeth, 8 diametral 
pitch, 1.75-in pitch radius, 20° pressure angle, 
m b =0.9. 




Fig. 14. Effect of rim compliance on crack 
propagation, m b = 1.3. (a) slotted gear, (b) partial- 
model gear. 
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Slotted 


n-0. 1 in 


u -0.0 1 in 



Fig. 15. Finite element model for study of rim/web 
compliance on crack propagation; 28 teeth, 

8 diametral pitch, 1.75-in pitch radius, 20° pressure 
angle. 


Table II. Tooth/rim deflections of uncracked gear (in), 
(Design: 28 teeth, 1. 75-in pitch radius, 8-diametral pitch) 


Backup 
ratio, mg 


r 

Slotted gear 

Web thickness 
u’=0.1 in u -0.01 in 

1.3 

Tooth 

0.000472 

0.000780 

0.005527 


Rim 

0.000279 

0.000506 

0.004494 

1.0 

Tooth 

0.000568 

0.000794 

0.005671 


Rim 

0.000375 

0.000536 

0.004732 

0.5 

Tooth 

0.001090 

0.000834 

0.005984 


Rim 

0.000849 

0.000598 

0.005134 


pitch, 1.75-in pitch radius, 20° pressure angle, and the 
standard fillet design. A plane stress, two-dimensional 
approximation was still used, but different thicknesses were 
specified for the tooth/rim face width, /, and the web 
thickness, w. Two different web thicknesses (w=0 . 1 and 
0.01 in) were studied and compared to the slotted baseline 
design. For all cases, the tooth face width was f= 0.25 in. As 
before, a tooth load of 2=364 lb was placed at the HPSTC, 
normal to the surface, and four hub nodes at the gear inner 
diameter were fixed to ground for boundary conditions. 
Table II gives the loaded tooth and rim deflections for an 
uncracked gear. Fig. 16 shows the effect of web thickness on 
crack propagation. Three backup ratios, m b = 1.3, 1.0, and 0.5, 
were studied. 

For m b =1.3, the vv=0.1-in model was slightly more 
compliant than the slotted model (Table II). However, tooth 
fractures occurred for all cases of initial crack location 
62 °< 00 ^ 119 ° for the w=0.1 -in model (Fig. 16b) while rim 
fractures occurred in the slotted model for 0^61° (Fig. 16a). 
This was also the trend for the m b = 1 .0 case except rim 
fractures occurred in the slotted model for %<83° (Fig. 16d), 
For both m*=1.3 and 1.0 and the vv=0.01-in model, the 
compliance was significantly increased (Table II) and rim 
fractures occurred for 6 0 <11° (Fig. 16c) and 0 0 <& 8° 
(Fig. 16f), respectively. For m b = 0.5, the slotted model and 
the vv=0.01-in model produced rim fractures for 0q<106 o 





Fig. 16. Effect of web thickness on crack propagation. 



Fig. 17. Gear tooth shapes, pitch radius, r =1.75 in. 
(a) 5.142857 pitch, 18 teeth, (b) 8 pitch, 28 teeth, 
(c) 16 pitch, 56 teeth. 


(Figs. 16g, 16i). The conclusion reached was that increased 
web compliance also lead to more rim fractures. However, 
when comparing a slotted gear to a webbed gear, compliance 
was not the only factor in determining tooth/rim fracture 
transition conditions. 

Effect of Gear Size on Crack Propagation 

The basic size of a tooth is determined by the 
fundamental equation: 



where P is the diametral pitch (in’ 1 ), N is the number of 
teeth, and r p is the pitch radius (in). Three different schemes 
were used in determining size effects on crack propagation: 
1) constant pitch radius, 2) constant number of teeth, and 
3) constant diametral pitch. 

Constant pitch radius. Fig. 17 shows three different 
tooth shapes for a constant pitch radius of r p =1.75in: 
a) 5.142857 diametral pitch, 18 teeth, b) 8 pitch, 28 teeth, 
and c) 16 pitch, 56 teeth. Case (b) was the baseline model 
described previously. Models for cases (a) and (c) are shown 
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Q 


r„= 1.75 


Fig. 18. Finite element models for gear size effect 
studies on crack propagation, (a) 5.142857 pitch, 
18 teeth, (b) 16 pitch, 56 teeth. 



Fig. 19. Effect of pitch and number of teeth on crack 
failure modes (constant pitch radius, r =1.75 in). 


in Fig. 18. Case (a) had 2226 elements and 6973 nodes while 
case (c) had 2396 elements and 7675 nodes. As before, the 
tooth load was placed at the HPSTC, normal to the surface, 
and four hub nodes at the gear inner diameter were fixed to 
ground for boundary conditions. Fig. 19 shows the design 
map (backup ratio and initial crack location effects on tooth 
or rim fractures) for the three cases. The 8-pitch, 28-tooth, 
model and the 16-pitch, 56-tooth model had nearly identical 
responses. The 5.142857-pitch, 18-tooth model had 
increased tooth-fracture conditions, indicating a decreased 
rim compliance condition. 

Constant number of teeth. Keeping the number of teeth 
constant and varying the size by proportionally varying the 
diametral pitch and pitch radius has the effect of scaling the 
design geometrically. That is, a design with N - 56 teeth, P=8 
(in' 1 ) diametral pitch, and ^=3.50 in is geometrically twice 
as big as a design with 77=56, P - 16 (in 1 ), and r p =1.75 in. 
Similarly, a A=56, P=5.283 (in 1 ), and r p =525 in is three 
times as big. To investigate this size effect on crack 
propagation, consider two simple machine elements: 
1) cantilever beam in bending, and 2) axially loaded bar 
(Pig. 20). A gear tooth can be roughly approximated by a 
cantilever beam in bending and an axially loaded bar in 



Fig, 20. Simple machine elements, (a) cantilever beam 
in bending, (b) axially loaded bar. 





Fig. 21. Gear tooth shapes, 8 pitch, (a) 28 teeth, 
r p = 1.-75 in, (b) 56 teeth, r p =3.50 in, (c) 84 teeth, 
r p -5 . 25 in. 

compression due to the tangential and radial tooth loads. The 
deflection, <5, and stress, <7, of a cantilever beam in bending 
with a rectangular cross section are: 

s _ 4Q X \ _ 6Qx ] 

° r 3 7 ° ~ 2 

tx 2*3 X 2 X 3 

where E is Young’s modulus. The deflection and stress of an 
axial bar in compression with a rectangular cross section are: 

8 <r = -2- 

Ex 2*3 X 2 X 3 

For both bending and compression, if the size is doubled, the 
magnitudes of the deflections are one-half of those for the 
original size for the same applied load and material as long 
as linear elastic conditions are applicable. Also, if the size is 
doubled, the magnitudes of the stresses are one-fourth of 
those for the original size. This proportioning is also 
applicable to the stress intensity factors, assuming linear 
elastic fracture mechanics. Recall that the crack propagation 
angles are a function of K n / K h Since Kj and K u are 
proportioned identically from the size effect, the ratio K n i K } 
remains the same. Therefore, the crack propagation paths are 
same. Thus, keeping a constant number of teeth and varying 
the diametral pitch and pitch radius has no effect on the 
crack propagation path. As a note, these results were- 
validated using finite element models and procedures 
described in the current study. 

Constant diametral pitch. Fig. 21 shows three different 
tooth shapes for a constant diametral pitch of P=8 in' 1 : a) 28 
teeth, />=1.75 in, b) 56 teeth, r p =3.50 in, and c) 84 teeth, 
r p =5. 25 in. Case (a) was the baseline model described 
previously. The model for case (b) was the same as in 
Fig. 18b considering the size effect of the previous section. 
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A model for case (c) was developed with 2448 elements and 
7969 nodes. Again, the tooth load was placed at the HPSTC, 
normal to the surface, and four hub nodes at the gear inner 
diameter were fixed to ground for boundary conditions. 
Fig. 22 shows the design map (backup ratio and initial crack 
location effects on tooth or rim fractures) for the three cases. 
The three cases had nearly the same response. Since the 
diametral pitch was identical, the basic tooth size (and thus 
rim size) was identical. Thus, the crack propagation paths 
were nearly the same. 

Effect of Pressure Angle on Crack Propagation 

A 20° pressure angle tooth and a 25° pressure angle 
tooth are depicted in Fig. 23. These were both for a 28-tooth, 
P - 8 diametral pitch, r p =1.75 in pitch radius design. Also, 
they had the same circular tooth thickness at the pitch point 
of k/ 2P=0.196 in. The 25° pressure angle tooth is wider at 
the base and narrower at the tip compared to the 20° pressure 
angle tooth. The design map for these cases is shown in 
Fig. 24. As seen, the 25° pressure angle tooth had increased 
tooth-fracture conditions, indicating a decreased rim 
compliance condition. 



Fig. 22. Effect of pitch radius and number of teeth on 
crack failure modes (constant 8-pitch gear). 




Fig. 23. Gear tooth shapes, 8 pitch, 28 teeth, pitch 
radius, r =1.75 in. (a) 20° pressure angle, 

(b) 25° pressure angle. 



Fig. 24. Effect of pressure angle on crack failure 
modes. 


Conclusions 

Design guidelines have been established to prevent 
catastrophic rim fracture failure modes when considering 
gear tooth bending fatigue. Analysis was performed using 
the finite element method with principles of linear elastic 
fracture mechanics. Crack propagation paths were predicted 
for a variety of gear tooth and rim configurations. The 
effects of rim and web thicknesses, initial crack locations, 
and gear tooth geometry factors such as diametral pitch, 
number of teeth, pitch radius, and tooth pressure angle were 
considered. The following conclusions were made: 

1) The occurrence of rim fractures significantly 
increased as the backup ratio (and thus, rim thickness) 
decreased. The occurrence of rim fractures also increased as 
the initial crack location was moved down the root of the 
tooth. A realistic design map of tooth/rim fracture modes 
included gear geometry, applied load, crack size, and 
material properties. 

2) Increased rim compliance increased the occurrence of 
rim fractures. Increased web compliance also increased the 
occurrence of rim fractures. When comparing slotted gears 
to web gears, however, compliance was not the only factor 
in determining tooth/rim fracture transition conditions. 

3) For gears with constant pitch radii, coarser-pitch 
teeth increased the occurrence of tooth fractures over rim 
fractures. Also, 25° pressure angle teeth increased the 
occurrence of tooth fractures over rim fractures when 
compared to 20° pressure angle teeth. For gears with 
constant number of teeth or gears with constant diametral 
pitch, varying size had little or no effect on crack 
propagation paths. 
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